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^ ■ Abstract 



C/3 



The recent measurement of a non-zero neutrino mixing angle #13 requires a modification of the 
tri-bimaximal mixing pattern that predicts a zero value for it. We propose a new neutrino mixing 
pattern based on a spontaneously-broken A4 flavor symmetry and a type-I seesaw mechanism. 
Our model allows for approximate tri-bimaximal mixing and non-zero #13, and contains a natural 



way to implement low and high energy CP violation in neutrino oscillations, and leptogenesis 
with a renormalizable Lagrangian. Both normal and inverted mass hierarchies are permitted 
within 3<r experimental bounds, with the prediction of small (large) deviations from maximality 
in the atmospheric mixing angle for the normal (inverted) case. Interestingly, we show that the 
inverted case is excluded by the global analysis in la experimental bounds, while the most recent 
MINOS data seem to favor the inverted case. Our model make predictions for the Dirac CP phase 
in the normal and inverted hierarchies, which can be tested in near-future neutrino oscillation 
experiments. Our model also predicts the effective mass |m ee | measurable in neutrinoless double 
beta decay to be in the range 0.04 < \m ee \ < 0.15 eV for the normal hierarchy and 0.06 < \m ee \ < 
0.11 eV for the inverted hierarchy, both of which are within the sensitivity of the next generation 
experiments. 
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I. INTRODUCTION 

The large values of the solar (612 — 35°) and atmospheric (6*23 — 45°) [l| neutrino mixing 
angles may be telling us about new symmetries in the lepton sector not present in the quark 
sector, and may provide us with a clue to the nature of the quark-lepton physics beyond the 
standard model. Theoretically, a great deal of effort has been put into constructing flavor 
models with high predictive power, especially those giving the tri-bimaximal (TBM) mixing 
angles 

013 = , 9 23 = j = 45° , 12 = sin" 1 ~ 35.3° . (1) 

However, the Daya Bay and RENO collaborations 3|,|4] have reported the first measurements 
of a non-zero value for the mixing angle #13: 

sin 2 20i3 = 0.092 ± 0.016(stat) ± 0.005(syst), (2) 

and 

sin 2 20i3 = 0.113 ± 0.013(stat) ±0.019(syst), (3) 

respectively, corresponding to an angle $13 ~ 9°. These results are in good agreement with 
the previous data from the T2K, MINOS and Double Chooz collaborations |5j. A non-zero 
value of 0i3 indicates that the TBM pattern for neutrino mixing should be modified. In ad- 
dition, at the Neutrino 2012 conference in Kyoto, the MINOS Collaboration has announced 
a non- maximal value for the atmospheric mixing angle 023 6J, 

sin 2 20 23 = 0.94^ ± 0.04 , (4) 

with maximal mixing disfavored at the 88% CL. This result, which was not included the 
global analysis in [3], comes from the analysis of disappearance in the MINOS accelerator 
beam, and points to one of two possible values for 2 3, namely 023 = 38° or 2 3 = 52°. If it 
holds, this result also calls for a deviation from the tri-bimaximal mixing pattern. 

Furthermore, the presence of CP violation in the lepton sector is still unknown. Experi- 
mentally, CP violation may become observable in a future generation of neutrino oscillation 



experiments [T2K, NOi^A] Theoretically, a flavor symmetry that describes and ex- 
plains the large reactor mixing angle 0i3 ~ 9° while keeping the TBM values 023 — 45° and 



9 12 — 35° may originate in two ways: (i) a large 6* 13 = Xc/y/2, with Ac the Cabbibo angle, 
mainly governed by higher-order corrections in the charged lepton sector 9j, where the TBM 
pattern is a good zero-order approximation to reality, or (ii) a large #13 from the neutrino 
sector itself through a new flavor symmetry without resorting to higher-order corrections in 



the charged lepton sector 10] . 



In this paper, we propose a new and simple model for the lepton sector with A 4 flavor 



symmetry in the framework of a type 



works using A4 flavor symmetries 
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seesaw mechanism. It is different from previous 



14f] 1 in that the Dirac neutrino Yukawa coupling 
constants do not all have the same magnitude. Our model can naturally explain the TBM 
large value of #13 and can also provide a possibility for low energy CP violation in neutrino 
oscillations with a renormalizable Lagrangian and small Yukawa coupling parameters, i.e. 
neutrino masses. The seesaw mechanism, besides explaining of smallness of the measured 
neutrino masses, has another appealing feature: generating the observed baryon asymmetry 



in our Universe by means of leptogenesis 



171 ] . Since the conventional A4 models realized 



with type-I or -III seesaw and a tree-level Lagrangian lead to an exact TBM and vanishing 
leptonic CP-asymmetries responsible for leptogenesis (due to the proportionality of the Y^Y V 
combination of the Dirac neutrino Yukawa matrix Y v to the unit matrix), authors usually 
introduce soft-breaking terms or higher-dimensional operators with many parameters, in 
order to explain the non-zero #13 as well as the non-vanishing CP-asymmetries. 

Our model is based on a renormalizable SU(2)l x U(l)y x A 4 Lagrangian with minimal 
Yukawa couplings, and gives rise to a non-degenerate Dirac neutrino Yukawa matrix and 
a unique CP-violation pattern. This opens the possibility of explaining the non-zero value 
of #i3 ~ 9° still maintaining TBM for the other two neutrino mixing angles #23 — 45° and 
9 12 — 35°; furthermore, this allows an economic way to achieve low energy CP violation in 
neutrino oscillations as well as high energy CP violation for leptogenesis. 

This paper is organized as follows. In the next section, we lay down the particle content 
and the field representations under the A4 flavor symmetry in our model, as well as explain 
the characteristic points of our model phenomenology at low and high energy. In Sec. Ill, we 
present the neutrino mixing angles, and how the low energy CP violation could be generated 



1 E.Ma and G.Rajasckaran [15| have introduced for the first time the A4 symmetry to avoid the mass 



degeneracy of \i and r under a fx-T symmetry [16( . 
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in both normal and inverted mass hierarchies, including our predictions for neutrinoless 
double beta decay. We give our conclusions in Sec. IV, and in Appendix A we outline the 
minimization of the scalar potential and the vacuum alignments. 



II. FLAVOR A 4 SYMMETRY FOR NON-ZERO 13 AND LEPTOGENESIS 

In the absence of flavor symmetries, particle masses and mixings are generally undeter- 
mined in a gauge theory. Here, to understand the present non-zero #13 and TBM angles 
(#12, #23) of the neutrino oscillation data and baryogenesis via leptogenesis, we propose a 
new discrete symmetry based on an A4 flavor symmetry for leptons in a renormalizable 
Lagrangian. 2 

The group A4 is the symmetry group of the tetrahedron, isomorphic to the finite group of 
the even permutations of four objects. The group A 4 has two generators, denoted S and T, 
satisfying the relations S 2 = T 3 = (ST) 3 = 1. In the three-dimensional real representation, 
S and T are given by 



'10 \ 



s 



\ 



^0 10^ 



T 



1 

1 



(5) 



0-10 

- 1 / V , 

A4 has four irreducible representations: one triplet 3 and three singlets 1, 1', 1". An A4 
triplet (di, 02, 03) transforms in the unitary representation by multiplication with the S and 
T matrices in Eq. (jSj) above, 
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a 2 



-a 2 
\~ a3 J 



T 



a i 
a 2 
\a 3 J 



(6) 



An A4 singlet a is invariant under the action of S (Sa = a), while the action of T produces 
Ta = a for 1, Ta = ua for 1', and Ta = u 2 a for 1", where u = e l2n ^ 3 is a complex cubic-root 
of unity. Products of two v4 4 representations decompose into irreducible representations 
according to the following multiplication rules: 3 <g> 3 = 3 S © 3 a © 1 © V © 1", V © 1" = 1, 



To include the quark sector, the symmetry could be promoted to the binary tetrahcdral group T' 18 1 
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1' <g) 1' = 1" and 1" ® 1" = 1'. Explicitly, if (ai, a 2 , 03) and (61, b 2 , 63) denote two triplets, 

(a®6) 3s = (a 2 6 3 + 0362,03^1 + ai^3,ai6 2 + a 2 &i) , 

(a g> 6) 3a = (a 2 b 3 - a 3 6 2 , a 3 &! - ai& 3 , - a 2 &i) , 

(a <g) 6)1 = aifei + a 2 b 2 + a 3 6 3 , 

(a ® 6)1/ = ai&i + ua 2 b 2 + w 2 ^^ , 

(a<E>6)i» = aifei + u 2 a 2 b 2 + wa 3 6 3 . (7) 

To make the presentation of our model physically more transparent, we define the T- 
flavor quantum number Tf through the eigenvalues of the operator T, for which T 3 = 1. In 
detail, we say that a field / has T-flavor Tf — 0, +1, or -1 when it is an eigenfield of the T 
operator with eigenvalue 1, u, u 2 , respectively (in short, with eigenvalue u Tf for T-flavor Tf, 
considering the cyclical properties of the cubic root of unity u). The T-flavor is an additive 
quantum number modulo 3. We also define the S'-flavor-parity through the eigenvalues of the 
operator S, which are +1 and -1 since S 2 = 1, and we speak of S'-flavor-even and S'-flavor- 
odd fields. For Aj-singlets, which are all S'-flavor-even, the 1 representation has no T-flavor 
(Tf = 0), the V representation has T-flavor Tf = +1, and the 1" representation has T-flavor 
Tf = —1. Since for ^U-triplets, the operators S and T do not commute, ^U-triplet fields 
cannot simultaneously have a definite T-flavor and a definite S'-flavor-parity. While the real 
representation of A4 in Eqs. (j3J), in which S is diagonal, is useful in writing the Lagrangian, 
the physical meaning of our model is more apparent in the T-flavor representation in which 
T is diagonal. This representation is obtained through the unitary transformation 



A A' = U u AUl, 



(8) 



where A is any A4 matrix in the real representation and 

(l 1 l\ 



1 id 2 u 
1 w u 2 



(9) 



We have 



/ 



S 



, 1 



-12 2 
2-12 
2 2-1 



T 



^10 0^ 

Ow 
co 2 



(10) 



Despite the physical advantages of the S' , T' representation, for clarity of exposition and 
to avoid confusion and complications, in this paper we use the real representation S, T 
almost exclusively. For reference, an A4 triplet field with components (a 4 , a 2 , 03) in the real 
representation can be expressed in terms of T-flavor eigenfields (a e , a M , a T ) (the notation 
comes from our lepton assignments below) as 



a e + a„ + a T a e + a; 2 a„ + ua T a e + ua^ + u 2 a T 

01 = V3 ' a2 = 71 ' a3 = 71 ' 



Inversely, 



_a 1 + a 2 + a 3 _ a x + oja 2 + ui 2 a 3 _ a x + tu 2 a 2 + ua 3 

We extend the standard model (SM) by the inclusion of an A 4 -triplet of right-handed 
S , f/(2)i-singlet Majorana neutrinos iV^, and the introduction of two types of scalar Higgs 
fields besides the usual SM S"?7(2)L-doublet Higgs bosons which we take to be an A4- 
singlet with no T-flavor (1 representation): a second S , [/(2)^-doublet of Higgs bosons 77, 
which is distinguished from $ by being an A 4 -triplet, and an 5'C/(2)i-singlet A 4 -triplet real 
scalar field x- 

•=(£)• %= (^)' »• j = w - (13) 

We assign each flavor of leptons to one of the three A4 singlet representations: the 
electron-flavor to the 1 (T-flavor 0), the muon flavor to the V (T-flavor +1), and the tau 
flavor to the 1" (T-flavor -1). (Note in this respect that our A 4 flavor group is not a 
symmetry under exchange of any two lepton flavors, like e and /1, for example. Our A 4 
flavor group is implemented as a global symmetry of the Lagrangian, later spontaneously 
broken, but some fields are not invariant under A 4 transformations, much in the same way 
as the implementation of SU(2)i x U(1)y in the SM, where left-handed and right-handed 
fermions are assigned to different representations of the gauge group.) Then we take the 
usual Higgs boson doublet $ to be invariant under A4, that is to be a flavor- singlet 1 with 
no T-flavor. The other Higgs doublet rj, the Higgs singlet x, and the singlet neutrinos Nr 
are assumed to be triplets under A4, and can so be used to introduce lepton-flavor violation 
in an A 4 symmetric Lagrangian. 

The field content of our model and the field assignments to SU(2)i x U(l)y x A 4 repre- 
sentations are summarized in Table HI These representation assignments and the requirement 
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that the Lagrangian be renormalizable and ^-symmetry forbid the presence of tree-level 
leptonic flavor- changing charged currents. 

The renormalizable Yukawa interactions in the neutrino and charged lepton sectors invari- 
ant under SU(2)l x U(1)y x A± are (including a Majorana mass term for the right-handed 
neutrinos) 

- C Yuk = yZL e {f}N R )i + y u 2 L^N R ) v + y^L T (fjN R ) v > 

+ y e L e <$> e R + y^L^ /i R + y T l T <$> r R + h.c, (14) 

where fj = IT2T1* and r 2 is a Pauli matrix. In this Lagrangian, each flavor of neutrinos and 
each flavor of charged leptons has its own independent Yukawa term, since they belong to 
different singlet representations 1,1', and 1" of A4: the neutrino Yukawa terms involve the 
A 4 -triplets r\ and N R , which combine into the appropriate singlet representation; the charged- 
lepton Yukawa terms involve the ^U-singlet $ and the ^U-singlet right-handed charged- 
leptons e R , fi R , and r R . The right-handed neutrinos have an additional Yukawa term that 
involves the ^-triplet SM-singlet Higgs \- The mass term ^M(N R N R )x for the right-handed 
neutrinos is necessary to implement the seesaw mechanism by making the right-handed 
neutrino mass parameter M large. 

The Higgs potential of our model contains many terms and is listed in Appendix A, 
Eqs. (IAip -( lA~7|) . 3 We spontaneously break the A4 flavor symmetry by giving non-zero vac- 
uum expectation values to some components of the ^-triplets x an d t\. As seen in Appendix 
A, the minimization of our scalar potential gives the following vacuum expectation values 



TABLE I: Representations of the fields under A 4 and SU(2)l X U(l)y- 



Field 


Lei L/j,, L T 


e R ,y R ,r R 


N R 


X 


$ 


77 




1, 1', 1" 


1, 1', 1" 


3 


3 


1 


3 


SU(2) L x U(1) Y 


(2,-|) 


(1,-1) 


(1,0) 


(1,0) 


(2,|) 


(2, J) 



3 We note that at TeV-scale the higher dimensional operators (d > 5) driven by x an d V fields are suppressed 
by a cutoff scale A which we assume is a very high energy scale, i.e. GUT or Planck scale. And in this 
paper we neglect the effects of higher dimensional operators. 
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(VEVs), all real: 

(^°> = ^0, (77?} = (r® = ( V ° 3 ) = ^= ? 0, ( X i)=v x ^O : (xa) = (X3> = 0. (15) 



The SM VEV v = (y/2G F y 1/2 = 246 GeV results from the combination v = Jv% + Zv*. 
The non-zero expectation value (</>°) = does not break the A A symmetry, because the 

standard model Higgs is ^-flavorless. The non-zero expectation value (77) = (v v , v v , v v )/y2 
breaks the S'-flavor-parity (771, 772, 773) — > (r]i, — 772, — ^3) but leaves the vacuum T-flavor Tf = 
0. In other words, after 77 acquires a non-zero VEV, the T-flavor is still conserved but the 
S- flavor-parity is not. Since 77 appears only in the Higgs sector and in interactions with the 
light leptons, we say that the light neutrino sector has a residual Z 3 symmetry expressed by 
the subgroup {1, T, T 2 } that leads to the conservation of T-flavor in terms involving mixing 
with the light neutrinos or interactions with the charged leptons. The non-zero expectation 
value (x) = (v x , 0, 0) maintains the S-flavor-parity of the vacuum (it is S'-flavor-even) but 
gives the vacuum the symmetric combination of T-flavors (ao + a + i +a_i) / \/3- That is, after 
X acquires a non-zero VEV, the S- flavor-parity is conserved but the T-flavor is not. Since x 
appears only in the Higgs sector and in interactions with the heavy Majorana neutrinos, we 
say that the heavy neutrino sector has a residual Z2 symmetry expressed by the subgroup 
{1, S} leading to the conservation of S'-flavor-parity in terms involving mixing or interactions 
with the heavy Majorana neutrinos. 

When a non-Abelian discrete symmetry like our v4 4 is considered, it is crucial to check 
the stability of the vacuum. In the presence of two y4 4 -triplet Higgs scalars \ and rj, Higgs 
potential terms involving both \ and 77, which would be written as V(xi]) in Eqs. ( 1A1I) - 
lJ]), would be problematic for vacuum stability. Such stability problems can be naturally 



solved, for instance, in the presence of extra dimensions or in supersymmetric dynamical 

cases, V(xr}) is not allowed or highly suppressed. 
The physical Higgs fields are obtained in the usual way. In the Higgs sector we have four 
Higgs doublets $, 771, r/ 2 an d 773, and three Higgs singlets Xii X2, an d Xs- They contain in 
total 16 degrees of freedom: six charged Higgs fields 77^23; with f]t = (v7)*i seven neutral 
Higgs scalars h, ^1,2,3, X1235 an< l three Higgs pseudoscalars ^1,2,3- We can write, after 
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electroweak- and /^-symmetry breaking and minimization of the potential, 

^ | > Xl = V X + Xl , X2 = X°2 , X3 = xl , 

^ I , J = 1,2,3. (16) 

^ (v„ + hj + iAj) 

The action of the residual Z 2 generator S on the physical fields is 

(N m , N m , N R3 ) -> (iV m , -iV ffi , -TVfla), (17) 

(xS,xS,xS)->(xS,-^,-x8), (is) 

(hi,h 2 ,h 3 ) ->■ (h 1 ,-h 2 ,-h 3 ), (19) 
(Ai,A 2 ,40->(Ai,-A 2 , -A 3 ), (20) 
^ ^ ) (^ -% + ; -% + ) ; (21) 

all other fields are invariant. The action of the residual Z 3 generator T on the physical fields 
is (the triplet fields ai, a 2 , and 03 and the triplet fields a e , a M , and a T are linear combinations 
of each other, see Eqs. (TITj)-(fT2l) 

(e,/i,r) -> (e,w/i,a; 2 r), (22) 
(u e , v^, v r ) -»■ (i/ e , (23) 
(iVife, iV^, iV^) -)• (JVife, wiVj^, w 2 iV fir ), (24) 
(xix%X°r)^(xico X %co 2 x% (25) 
(/i e , h T ) ->■ (/i e , w 2 /i r ), (26) 
(A e , A M , A T ) -»• (A e , wij,, w 2 A r ), (27) 

(^ Vi, Vr) Wr £> W V )> (28) 

all other fields are invariant. 

After electroweak and A 4 symmetry breaking, the neutral Higgs fields acquire vacuum 
expectation values and give masses to the charged-leptons and neutrinos: the Higgs doublet 
gives Dirac masses to the charge leptons, the Higgs doublet 77 gives Dirac masses to the three 
SM neutrinos, and the Higgs singlet x gives a Majorana mass to the right-handed neutrino 

Nr. 
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The charged lepton mass matrix is automatically diagonal due to the /U-singlet nature 
of the charged lepton and SM-Higgs fields. The right-handed neutrino mass has the (large) 
Majorana mass contribution M and a contribution induced by the electroweak-singlet A4- 
triplet Higgs boson \ when the ^-symmetry is spontaneously broken. 

After the breaking of the flavor and electroweak symmetries, with the VEV alignments 
as in Eq. (fl5|) . the charged lepton, Dirac neutrino and right-handed neutrino mass terms 
from the Lagrangian fTHj) result in 



1 v ( 

- C m = ^= (y e e L e R + y^ L fi R + y T f L r R ) + -|| {y\v e + y^v^ + yp T ) N R1 



+ [y u x v e + y 2 uv^ + V) A^ 2 + u e + y 2 u 2 v^ + j/£wP T ) iV ra 



^ ( atc at_ 1 ?\rc 7\r 1 Arc j\r n 1 yR v x, 



+ y + + N R3 N m ) + ^(JV&JVa, + N R3 N R2 ) + h.c. . (29) 

This form shows clearly that the terms in v v break the S'-flavor-parity symmetry (|T7j) - (l2Tl) . 
while the other mass terms preserve it. Passing to the T-flavor eigenfields 

N RC = ^ + + (30) 
, r N R1 + uN R2 + u 2 N R3 

^Rfi = , (31J 

^Rt ) l^J 

with respective T-flavor Tf = 0, +1, —1, the lepton mass Lagrangian reads 

- C m = -j= {y e eLe R + V/MpLfiR + y T f L T R ) 
\/2 

+ (ViVeNRe + J/^V^ + yp T N Rr ) 

+ fiWe^Re + W^Rr + Wr^R,) 

Vr v x 



N Re N Re + N R ^N R , + N c Rt N Rt 

~l(We+W,+Wr) (NRe + Nr, + N Rt ) ] + h.C. . (33) 

This form shows clearly that the terms in v x break the T-flavor symmetry (j22p -(j28 p . while 
the other mass terms preserve it. 

Inspection of the mass terms in Eq. (1331) indicates that, with the VEV alignments in 
Eq. ()15p . the A 4 symmetry is spontaneously broken to a residual Z 2 symmetry in the heavy 
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Majorana neutrino sector (conservation of S'-flavor-parity in terms not involving v v or ^1,2,3) 
and a residual Z3 symmetry in the Dirac neutrino sector (conservation of T-flavor in terms 
not involving v x or Xi)- 

The mass terms in Eq. ( 129"j) and the charged gauge interactions in the weak eigenstate 
basis can be written in (block) matrix form as, using N R mi)V c L = VZm^Nn, 

- C m w = hr R M R N R + Wm D N R + T L m t i R + ^W;T L ^v L + h.c. (34) 
= \ (- L W R \ I ° ^ J f Vl I + W« + -^W^T^z + h.c. (35) 
Here £ = (e, /1, r), z/ = (u e , u T ), N R = (N m , N R2 , N m ), and 



7! 



V 0^ 

y M 

^0 y r y 



(36) 



m - Vv Y - - v " 



( 



M, 



R 



\ 



(37) 



MO 

M Vr v x ■ (38) 
V° V&X M I 

To find the neutrino masses and mixing matrix we are to diagonalize the 6x6 matrix 



(39) 



We start by diagonalizing M R . For this purpose, we perform a basis rotation N R = U R N R , 
so that the right-handed Majorana mass matrix M R becomes a diagonal matrix M R with 
real and positive mass eigenvalues Mi = aM, M2 = M and M3 = bM, 




I 



M R = U R M R U R = MU T R 



1 

1 Ke^ 
«e* 1 



aM 
M 
bM 



(40) 
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where k = \y R v x /M\ and £ = &rg(y R v x /M). We find a = \/l + k 2 + 2k cos £, 6 
-y/l + — 2k cos£, and a diagonalizing matrix 

\ 



( V2 V e'T ( 



V2 



1 -1 
1 1 



1 




(41) 



with phases 



tpi = tan -1 ( ) and tp2 = tan _] 



(42) 



. 1 + K COS £ / VI — K COS £ / 

As the magnitude of k defined in Eq. ( )40|) decreases, the phases ipxp go to or ir. At this 
point, 



-C 



mW 



W N' 



R 



f e L m t e R + -^W-lL-fvL + h.c. (43) 
V2 



m T D M R I \N R 



With TYl£) = TTIdUr. 

Now we take the limit of large M (seesaw mechanism) and focus on the mass matrix of 
the light neutrinos M u , 

-£mw = )^7ZM v v c L + hmdR + ■A=W~I£'fu L + h.c. + terms in N R (44) 

with 



M u = -m D M R l m T D . 



(45) 



We perform basis rotations from weak to mass eigenstates in the leptonic sector, 



h = p;e L , £ R = p;i r , v L = ulpy L , 



(46) 



where P% and P v are phase matrices and U„ is a unitary matrix chosen so as the matrix 



m 



„ = UlPtMyPfUt = -UlPtm D U R M^{T]lPtm D U R ) 



(47) 



is diagonal. Then from the charged current term in Eq. (|43p we obtain the lepton mixing 
matrix C/pmns as 



t^PMNS — PpPyUy. 



(48) 
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The matrix £7pmns can be written in terms of three mixing angles and three CP-odd phases 
(one for the Dirac neutrinos and two for the Majorana neutrinos) as [l| 

/ ,. ,. ,. .. . ,. , '•><•;• \ 



PMNS 



C13C12 C13S12 si 3 e l5c 

~C23Si2 — •523Ci2Si3e t(Scp C23C 12 — S23Si2Si3e iScp S23 c 13 



S23 s 12 — C23Ci2Si3e 



~ s 23 c 12 — C23S 12 Si3e 



"CP 



C23 c 13 



Qu, (49) 



/ 



where Q v = Diag(e~ 4ipi / 2 , e _ ^ 2//2 , 1), and = sin#y and Cy = cos 6^-. 

It is important to notice that the phase matrix P u can be rotated away by choosing 
the matrix Pu = P v , i.e. by an appropriate redefinition of the left-handed charged lepton 
fields, which is always possible. This is an important point because the phase matrix P v 
accompanies the Dirac- neutrino mass matrix rhp and ultimately the neutrino Yukawa matrix 
Y v in Eq. (j3"7|) . This means that complex phases in Y v can always be rotated away by 
appropriately choosing the phases of left-handed charged lepton fields. Hence without loss 
of generality the eigenvalues y\, y%, and y% of Y v can be real and positive. The Yukawa 
matrix Y v can then be written as 



^ 2/i 0^ 



\ 



y 2 
1 



(50) 



where y x = \y\lvl\,y2 = IvX/Usl and U u is given in Eq. ©. 

Concerning CP violation, we notice that the CP phases ipi,ip2 coming from Mr only 
take part in low-energy CP violation, as can be seen in Eqs. P0H50|) . Any CP- violation 
relevant for leptogenesis is associated with the neutrino Yukawa matrix Y v = Y v Ur and the 
combination of Dirac neutrino Yukawa matrices, H = YjY u = U^YJY^Ur, which is 



/ 



h = m 



v\2 



l+%f+y| 



V 



V2 



-vl 



■ ^21 

m^{y\ 



3 e^(i/2 



1) 



i) l(i + vi) 



(51) 



/ 



where ipij = ipi—ipj. As expected, in the limit 



M\ 



\y% 1 , i.e. 2/1,2 —> 1, the off-diagonal 



entries of H vanish, and there is no CP violation useful for leptogenesis. If the Dirac neutrino 
Yukawa couplings y±, y%, and y\ differ in magnitude, they can play a role in baryogenesis via 
leptogenesis and non-zero #13 ~ 9° with TBM (#23 — 45°, # 12 — 35°). Therefore, a low energy 
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CP violation in neutrino oscillation and/or a high energy CP violation in leptogenesis can 
be generated by the non- degeneracy of the Dirac neutrino Yukawa couplings and a non-zero 
phase £ coming from Mr. 

In the following section we investigate the low energy phenomenology, namely the pos- 
sible values of the light neutrino mixing angles, how the low energy CP violation could be 
generated in both normal and inverted mass hierarchies, and neutrinoless double beta decay, 
which is a probe of lepton number violation at low energy. 



III. PHENOMENOLOGY OF LIGHT NEUTRINOS 



After seesawing, in a basis where charged lepton and heavy neutrino masses are real and 
diagonal, the light neutrino mass matrix is given by 

„,2 



m. 



( (i + 2 -^)vl (1-^)^2 



m 



m 



\ 



)ym (i 



)y\ 



\ 



'1 e'*l\ 



e'^l _ 3e"^2 
2a 2b 
e^i , 3e'^2 s 

-2a~ + -^rm 



'1 -i- e 

- L ~ l ~ 2a 



2I> 



J 



(52) 



where = v v Y v / a/2 and we have defined an overall scale m = v^\y^ 2 \/(6M) for the light 
neutrino masses. The mass matrix m v is diagonalized by the PMNS mixing matrix £/pmns 
as described above, 



^pmns Diag(mi, m 2 ,m 3 ) Up MNS . 



(53) 



Here rrii (i — 1, 2, 3) are the light neutrino masses. As is well-known, because of the observed 



hierarchy \Am\ tr 



m 



-m\\ Amg ol = 777-2— m i > 0) an d the requirement of a Mikheyev- 



Smirnov-Wolfenstein resonance for solar neutrinos, there are two possible neutrino mass 
spectra: (i) the normal mass hierarchy (NMH) mi < m 2 < m 3 , and (ii) the inverted mass 
hierarchy (IMH) 7773 < mi < 7772. 

Interestingly, the combination U^Ur in Eq. fl52l) reflects an exact TBM: 

( [2 \ n \ / j& n q \ 



3 v 7 ^ 

1 1 

1_ J_ 1 



1 

>/2 



1 



e 



(54) 
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Therefore Eq. (|52|) directly indicates that there could be deviations from the exact TBM 
if the Dirac neutrino Yukawa couplings do not have the same magnitude. In the limit 
1 2/2! = l?/3l (2/2 — ► 1); the mass matrix in Eq. ( 15 2 p acquires a /x-r symmetry that leads to 
# 13 = and 6*23 = — vr/4. Moreover, in the limit \y"\ = \y 2 \ = \y%\ (yi,yz — > 1), the mass 
matrix (I52p gives the TBM angles in Eq. ([I]) and the corresponding mass eigenvalues 

3m 3mo 



mi 



m 2 = 3mo 



m 3 



(55) 



a b 
These mass eigenvalues are disconnected from the mixing angles. However, recent neutrino 
data, i.e. # 13 7^ 0, require deviations of y^ 2 from unity, leading to a possibility to search 
for CP violation in neutrino oscillation experiments. These deviations generate relations 
between mixing angles and mass eigenvalues. 

To diagonalize the above mass matrix Eq. (|52|) . we consider the hermitian matrix m„mj = 
^pmns Diag(m^, m 2 , 777.3) ^pmns> fr° m which we obtain the masses and mixing angles. To see 
how the neutrino mass matrix given by Eq. (l52l can lead to deviations of neutrino mixing 
angles from their TBM values, we first introduce three small quantities (i = 1,2,3), 
which are responsible for the deviations of the Ojk from their TBM values: 

7T / 1 

#23 = - ^ + d , 013 = £2 , O12 = sin" 1 I — 
Then the PMNS mixing matrix up to order e, can be written as 

v. 

f^PMNS = 



£3 • 



(56) 



V 



-l+ei+e 3 V2 



e 2 e l0 CP 



V3 



V3 

V6 



£2e 



J^-o^CP 

V6 6 



1+ei 



Q„ + Oft) . (57) 



The small deviation ei from the maximality of the atmospheric mixing angle #23 is ex- 
pressed in terms of the parameters in Eq. (IBip in Appendix B as 

R(l + y 2 )-S(y 2 -l) 



tanei — , . 

i?(y 2 - 1) - 5(1 + jfe) 

In the limit of 2/1,2/2 ~ ^ 1> e i g° es to zero (or equivalently #23 _ > due to R, S 

The reactor angle 6*13 and the Dirac-CP phase 5c p are expressed as 

2/i|s 23 (3Q - P)y 2 - c 23 (3Q + P) - 3i{s 23 (R - S)y 2 + c 23 (R + 5)}| 



(58) 
->■ 0. 



tan26> 13 
tan^p 



(F + G+f + f)(c| 3 



?/2 S 23) 



j, 2 (F + G-5f) sin 2^3 



(P + Q)(fl-S)-j/2(P-Q)0R+S) 



(59) 
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where the parameters P,Q,F, G, K, D and A are given in Eq. f ]Blj) in Appendix B. In the 
limit of jji, y2 —> 1, the parameters Q, R, S go to zero, which in turn leads to 6* 13 — > and 
5c p —> as expected. Finally, the solar mixing angle is given by 

2/ 2 c 23 (3Q - P) + s 23 (3Q + P) 



tan26> 



12 



2yi 



(60) 



Ci 3 (*2-*l) 

Since in the limit 2/1,2/2 — > 1 the parameters in Eq. (1601 behave as Q — > 0,P — > 6(s — 
1), — > 3(1 + t) and ^2 — > 6(1 + ^2), it is clear that the mixing angle tan26*i2 goes to 
2V2, that is, 0i2 -> Bin -1 (l/V3). 

The squared-mass eigenvalues of the three light neutrinos result in 







m\ 


= mgj 






ml 


= mjj| 






m 2 3 


= ™o{ 







S? 2 * i + - Vl 



2/2C 23 (3Q-P) + S23(3Q + P) 



2c 



sin 26* 12 
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2 xT , 2 xT , ^ y 2 c 23 (3g-p) + g23 (3g + p) . 

m <( c 12 *i + s 12 tf 2 + 2/i r sm 26>i 2 

^Ci 3 

F + G + ^ + ^) (4 + 2/ 2 2 4) + 2/2 (p + G ~ ^) Bin 2ft 



23 



-13 



2 _ 2/1 sin26>i 3 

13 o 



c 23 ((3Q + P) cos 5 C P - 3{R + S) sin S CP ) 



S23IJ2 ((30; - P) cos 5cp + 3(P - S) sin 5 CP ) | . 



(61) 



We see from Eqs. ( 16"U|) and (|6ip that the deviation e 3 from tri-maximality of solar mixing 
angle #i 2 can be expressed as 

32/im2{2/2C 23 (3g - P) + s 23 (3Q + P)} 



2-\/2 cos 2e 3 + sin 2e 3 



Ci 3 Am 2 



(62) 
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Now we perform a numerical analysis using the linear algebra tools in Ref. [20j. The Daya 
Bay and RENO experiments have accomplished the measurement of three mixing angles 
#i2,#2 3 , and #13 from three kinds of neutrino oscillation experiments. A combined analysis 
of the data from the T2K, MINOS, Double Chooz, Daya Bay and RENO experiments 
shows [7] that, for the normal mass hierarchy (NMH) and inverted mass hierarchy (IMH), 
respectively, 



tin' ft, = 0.026«;S£ <«SS NMH, 



or equivalently 



'13 



9.28°! °S: tiZ] NMH, 



n n97 +0 - 003 (+°- 01 °) TMW 
0-027 -0.004 (-0.011) 1MH 



9.46-;- [SSj IMH 



(63) 



(64) 
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at the la (3a) level. The hypothesis #13 = is now rejected at the 8a significance level. In 
addition to the measurement of the mixing angle 6*13, the global fit of the neutrino mixing 
angles and of the mass-squared differences at the la (3a) level is given by Q| 



a _ oa ^+0.92° (+3.02°) ft _ 4q o+4.60° (+8.70°) ivrA/rrr 

V 12 - 34.45 _ 105o (_ 3 _ 14 o) , #23 - 44.43 _ 28JO ( _ 5 7go) IN MM, 



AR 70 o+2.89° (+6.41°) tt\ /Ttt 

46.72 _ 401o ( _ 8 07o) IMH 



? .0+0.08 (+0.24) NMR 

Am 2 [10" 5 eV 2 l - 7 62 +0 - 19 (+a58) Am 2 [10" 3 eV 2 l - I ~o.io (-0.27) > mvin , fi ^ 
Am Sol [lU eV j - f.b2_ 019 ( _ 050) , Am Atm [lU eV j - < +010 (+028) {^) 

y 2 - 4U -o.o7 (-o.25) ' iivm 

The matrices mjj and Mr in Eq. (|52|) contain seven parameters : y%, M, v^, yi, y 2 , k, £. 
The first three (7/3, M, and v,,) lead to the overall neutrino scale parameter rriQ. The next four 
(2/1 5 2/2 5 give rise to the deviations from TBM as well as the CP phases and corrections 
to the mass eigenvalues (see Eq. (1551) ). 

In our numerical examples, we take M = 10 TeV and v v = v $ = 123 GeV, for simplicity, 
as inputs. Since the neutrino masses are sensitive to the combination m = t> 2 |?/3 2 |/(6M), 
other choices of M and v v give identical results. Then the parameters mo, Z/i, 2/2; £ can be 
determined from the experimental results of three mixing angles, #12, #13, #23, and the two 
mass squared differences, Am^, Am^. In addition, the CP phases 8c p, y?i,2 can be predicted 
after determining the model parameters. Using the formulae for the neutrino mixing angles 
and masses and our values of M, v v ,v$, we obtain the following allowed regions of the 
unknown model parameters: for the normal mass hierarchy (NMH) 4 , 

0.17<k<0.90, 0.74 < Vl < 1.0 , 0.90 < u 2 < 1.11 , 

94 < £ < H9° 

~ , 1.8<m xl0- 2 [eV]<6.0; (66) 

240° < i < 265° 

for the inverted mass hierarchy (IMH), 

0.31 < k < 0.92 , 0.84 < Vl < 1.15 , 0.65 < y 2 < 1.28 , 

1.7<m xlO- 2 [eV]<4.5. (67) 



90°<e< H7° 



245° < i < 265° 

Note that here we have used the 3a experimental bounds on # 12 , #23, Am 2 1; Am^ in Eq. (1651) . 



except for #13 < 12° for which we use the values in Eqs. (I66f67l) . For these parameter regions, 



4 When |/2 = 1 and around there, there exist other parameter spaces giving very small values of #13. So, 
we have neglected them in our numerical result for normal mass hierarchy. 
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FIG. 1: The reactor mixing angle #13 versus the ratio of first-to-second generation neutrino Yukawa 
couplings 2/1/2/2 (teft plot) and the parameter k = \y u R v x /M\ (right plot). The (red) crosses and 
(blue) dots represent the results for the normal and the inverted mass hierarchy, respectively. The 
horizontal solid (dotted) lines in both plots indicate the upper and lower bounds on #13 for inverted 
(normal) mass hierarchy given in Eq. (|64p at the 3a level. 




FIG. 2: The atmospheric mixing angle #23 versus the phase £ of the parameter combination 
y R v x /M. The (red) crosses and (blue) dots represent the results for the normal and inverted 
mass hierarchy, respectively. 

we investigate how a non-zero #13 can be determined for the normal and inverted mass 
hierarchy. In Figs. HUH the data points represented by blue dots and red crosses indicate 
results for the inverted and normal mass hierarchy, respectively. The left-hand-side plot 
in Fig. CD shows how the mixing angle 6*13 depends on the ratio yx/y2 = Vx/Vi °f the first- 
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FIG. 3: The Jarlskog invariant Jcp versus the reactor angle #13 (left plot), and the Dirac CP phase 
5c p versus #23 (right plot). The (red) crosses and (blue) dots represent the results for the normal 
and inverted mass hierarchy, respectively. The vertical solid (dashed) lines in both plots indicate 
the upper and lower bounds on #13 for the inverted (normal) mass hierarchy given in Eq. (|64p at 
the 3a level. 

and second-generation neutrino Yukawa couplings; the right-hand-side plot shows how 13 
depends on the parameter k = \y^v x /M\. We see that the measured value of #13 from the 
Daya Bay and RENO experiments can be achieved at 3<r's for 0.75 < yi/y2 < 1 (NMH), 
1.1 < < 1.3 and y x /y 2 ~ 0.9 (IMH), 0.17 < k < 0.82 (NMH) and 0.3 < k < 0.74 

(IMH). Fig. HJshows the atmospheric mixing angle 623 as a function of the phase £ of y u R v x /M. 



To see how the parameters are correlated with low-energy CP violation observables mea- 
surable through neutrino oscillations, we consider the leptonic CP violation parameter de- 



fined by the Jarlskog invariant 



2l| 



1 



J C p = ImlUeiU^U^U*^ = - sin 26 12 sin 26> 23 sin 26 13 cos 6 13 sin 5 C p- 



(68) 



The Jarlskog invariant Jcp can be expressed in terms of the elements of the matrix h 

Im{/ii 2 /i 2 3^3i} 



CP 



(69) 



The behavior of Jcp as a function of #13 is plotted on the left plot of Fig. |3j We see that 
the value of \ J C p\ lies in the range - 0.04 (NMH) and 0.02 - 0.04 (IMH) for the measured 
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value of #13 at 3<t's. Also, in our model we have 



Im{/ii 2 /l23^3l} 



27m£ 



ViVii 1 -2/2) sin ^ 2 {--} , 



(70) 



4a 4 6 3 



in which { } stands for a complicated lengthy function of yi, 2/2, cl, b, ip\ and ip2- Clearly, 

Eq. flTOjl indicates that in the limit of y 2 — > 1 or simp 2 the leptonic CP violation Jqp 
goes to zero. When 1/2 7^ 1, he. for the normal hierarchy case, Jcp could go to zero as sin ip2 
of Eq. ( 17D1) . In the case of the inverted hierarchy, Jqp has nonzero values for the measured 
range of #13 while Jcp goes to zero for # 13 — > 0, which corresponds to y 2 — > 1. The right plot 
of Fig. [3] shows the behavior of the Dirac CP phase 5c p as a function of 823, where Sep can 
have discrete values around 50°, 120°, 230° and 310° for the inverted mass hierarchy (for the 
normal mass hierarchy, 5c p can vary over a wide range except near 90° and 270°). Future 
precise measurements of #23, whether #23 > 45° or #23 < 45°, will provide more information 
on 5cp- 

Fig. H] shows how the values of #13 depend on the mixing angles 823 and 812 . As can be seen 
in the left plot of Fig. HI the behavior of 6*23 in terms of the measured values of #13 at 3cr's for 
the normal hierarchy is different than for the inverted hierarchy. For the normal hierarchy 
we see that the measured values of 8i 3 can be achieved for 43° < 823 < 47° and 6*23 7^ 45°, 
with small deviations from maximality, while for the inverted hierarchy 50° < 823 < 53.1° 
and 38.6° < 823 < 40°, which are excluded at la by the experimental bounds as can be 
seen in Eq. ( 165]) . 5 From the right plot of Fig. HI we see that the predictions for 8 13 do not 
strongly depend on # 12 in the allowed region. 

Moreover, we can straightforwardly obtain the effective neutrino mass \m ee \ that charac- 
terizes the amplitude for neutrinoless double beta decay : 



where £/pmns is given in Eq. (1571) . The left and right plots in Fig. [5] show the behavior of the 
effective neutrino mass \m ee \ in terms of 6^3 and the lightest neutrino mass, respectively. In 
the left plot of Fig. [5], for the measured values of 813 at 3er's, the effective neutrino mass \m ee \ 



can be in the range 0.04 < |m ee |[eV] < 0.15 (NMH) or 0.06 < |m ee |[eV] < 0.11 (IMH). The 



5 Interestingly, the most recent data of MINOS seem to disfavor the maximal mixing in the atmospheric 
mixing angle, Eq. (j4]), indicating that the inverted mass hierarchy may be favored. 




(71) 
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FIG. 4: The behaviors of #23 and O12 in terms of #13. The red crosses and the blue dots represent 
results for the normal mass hierarchy and the inverted mass hierarchy, respectively. The solid 
(dashed) vertical lines represent the experimental bounds of Eq. (|65p at 3cr's for the inverted 
(normal) mass hierarchy. The horizontal dotted lines indicate the la experimental bounds in 
Eq. (J65])- 




8 13 [Deg.] m llgh , esl [eV] 



FIG. 5: Plots of \m ee \ as a function of #13 and mutest- The red crosses and the blue dots represent 
results for the normal and the inverted mass hierarchy, respectively. The vertical solid (dashed) 
lines show the experimental bounds of Eq. (|65p at 3<r's for the inverted (normal) mass hierarchy. 

right plot of Fig. |5] shows | function of raughtest, where rniightest = vtt-i for the normal 

mass hierarchy and mi ig htest — ^3 for the inverted mass hierarchy. Our model predicts that 
the effective mass \m ee \ is within the sensitivity of planned neutrinoless double-beta decay 
experiments. 
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IV. CONCLUSIONS 



We have suggested a novel and simple scenario to generate neutrino masses and mixings 
with a discrete A 4 symmetry that is spontaneously broken. In particular our model can ac- 
commodate in a renormalizable Lagrangian a large value of the mixing angle, #13, consistent 
with the recent reactor neutrino experiments Daya Bay and RENO, as well as high energy 
CP violation interesting for leptogenesis. 

In our model we have introduced a right-handed neutrino Nr, a real gauge-singlet scalar 
X, and an SU (2) ^-doublet scalar 77, all of which are A A triplets. The light neutrino masses 
are generated by a seesaw mechanism in which we have assumed the right-handed neutrino 
masses are at the TeV scale (to evade the introduction of higher dimensional operators). 
Getting VEVs along the direction (%} = v x (l,0, 0) and (77°) = v v (l, 1, 1), which break the 
A4 symmetry down to a Z 2 (S- flavor-parity) and a Z 3 (T-flavor) symmetry, respectively, 
one obtains bimaximal mixing at the right-handed neutrino sector and trimaximal mixing 
at the light Dirac neutrino sector with non-degenerate Yukawa couplings that deform the 
exact TBM pattern. The resulting light neutrino mixing matrix is in the form of a deviated 
TBM generated through unequal neutrino Yukawa couplings, as can be seen in Figure [TJ 
In the limiting case of equal active-neutrino Yukawa couplings, the mixing matrix recovers 
the exact TBM. In addition, we have shown that unequal neutrino Yukawa couplings can 
provide a source of high-energy CP violation, perhaps strong enough to be responsible for 
leptogenesis. The stability of the vacuum alignments we assume are guaranteed, for example, 
by embedding our model in an extra dimension. 

We showed that deviations from the TBM of about 20% are enough to explain # 13 ~ 9°. 
We predicted that the CP violating Dirac phase Sep may have discrete values (see Figure [3]). 
Therefore the measurement of the phase Sep in the next generation neutrino experiments 
can rule out or support our model. We have also shown that the inverted mass hierarchy 
may be excluded by a global analysis using la experimental bounds, while the most recent 
MINOS data seem to favor it. We also predicted an effective neutrino mass in neutrinoless 
double-beta decay in the range, 0.04 < |m ee |[eV] < 0.15 (for the normal hierarchy) and 
0.06 < |m ee |[eV] < 0.11 (for the inverted hierarchy), both ranges within reach of near-future 
neutrinoless double beta decay experiments. 
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Appendix A: The Higgs potential 



In this Appendix we present our Higgs potential and its minimization, as well as our 
prescription for effecting the stability of the vacuum alignment. We solve the vacuum align- 



ment problem by extending the model into a spatial extra dimension y [12|. We assume that 
each field lives on a 4D brane either at y = or at y = L, as shown in Fig. |6j The heavy 
neutrino masses arise from local operators at y — 0, while the charged fermion masses and 
the neutrino Yukawa interactions are realized by non-local effects involving both branes A 
rigorous explanation of this possibility is beyond the scope of this paper. 

The most general renormalizable scalar potential for the Higgs fields <3>, rj and x, invariant 
under SU{2) L x U(l)y x A± and obeying the conditions in the previous paragraph, is then 
given by 



V = V y=0 + V y=L , (Al) 



where 



V y=0 = V($) + V(r]) + V(rj$) , (A2) 
V y=L = V(x) , (A3) 



and 



V(rj) = fi 2 v {v j vh + Xltfvhtfvh + ^v)i'Wv)i" + *Uv t v)s.(tfv)s. (A4) 
V{$) = //|($t$) + A*($t$)2 ; 

V(x) = f4(xx)i + X i(xx)i(xx)i + ^2(xx)v(xx)v' + >$(xxhs(xx)3 a (A5) 

+ ^4(xx)sa(xx)sa + ^Uxx)s s (xx)s a + £iX(xx)s. + &x(xx)s a , 
V{r)$) = \f{7i*7])i{&$) + Af [(7 7 t$)($t 77 )] 1 + | A f [(^t$)( ?7 t$)] 1 + h . c X ( A6 ) 

+ {Af (Jv)*.(J*) +h-c.} + {Af (i^W*) + h.c.} . (A7) 

Here fi v , and have mass dimension-1, while A^ 5 , A*, A* 5 and Af 5 are 

dimensionless. In V(r/$) the usual mixing term &rj is forbidden by the A4 symmetry. In 
the scalar potential flAIj) -( lA7|) we have for simplicity assumed that CP is conserved, and 
the couplings Af , Af and A 5 are real. 
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FIG. 6: Fifth dimension y and locations of scalar and fermion fields on the brane at y = and 
y = L. 

The vacuum configuration is obtained by the vanishing of the derivative of V with respect 
to each component of the scalar fields $>,r]i,Xi (i = 1,2,3). The vacuum alignment of the 
field rj is determined by 



V2 



dV 



V2 



V2 



dr]i (v°) 



dV 



dv2 (»»?>= 
dV 



= V -f{2vl(\l + Ag) + 2nl + « + <)(2A? - \ v 2 + 4A§) 

+ 4(Af + Af + 2Af )} + 3v V2 v m v*\f = , 

= ^{2<(A? + A3) + 2^ + « + <)(2A? - XI + 4Ag) 

+ 4(Af + Af + 2Af )} + 3v Vl v m v*\f = , 

= V ~f{K( X i + X D + 2/4 + « + <)(2A? + 4Ag) 
+ 4(Af + Af + 2Af )} + 3v m v m v^\f = . 
From this set of three equations, we obtain the solution 

fa?> = <%°> = <%°> = ^ 

-3v*Af ± x /94Af 2 - 2(3A? + 4A|DM + v |(Af + Af + 2Af )) 



^3 



(A8) 



2(3A? + 4A 3 ? ) 
This VEV breaks A± down to a residual Z3. 

The vanishing of the derivative of V with respect to $ reads 



^ . (A9) 



y/2 



5$o 



{»|A* + 4 + i(Af + Af + 2Af )« + ^ + «J,)} 



+ 3 1^ t>,j 2 A 4 



)/<& 



. 



(A10) 
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The real- valued solution of Eq. ( lAlOj) . for real- valued parameters, is 



9a 2 c+ W3(45 3 6 3 + 27a 4 c 2 ; 



1/3 



9a 2 c + x /3(4a 3 6 3 + 27a 4 c 2 ; 

^(is) 1 ^ 



1/3 



(All) 



where a = A*, 6 = /x| + |vJ(Af + Af + 2Af ) and £ = 3w 3 Af . 

Finally, the minimization equations for the vacuum configuration of x are given by 

dV 



dxi 
dV 



dx2 
dV 



9x3 



(Xi)=v Xi 



2^ X1 (ft + (2A* - A* + 4A£)(t£, + t&) + 2(A? + AJ)^) 
6^X2^X3 = , 

2v X2 (fi + (2Af - A* + 4Af)(^ + t&) + 2(A? + A*)t&) 



+ G$v xi v xs = 



2 %3 + (2A? - A* + 4A*)« + vl) + 2(A* + Aj)i£.) 



(A12) 



From these equations, we obtain the solution 6 



(Xi) 



2(Af + a*; 



7^ , ( X2 ) = (X3> = . 



(A13) 



Appendix B: Parametrization of the neutrino mass matrix 



We parametrize the hermitian matrix m u m\, as follows: 



/ 



m v m\ = r«Q 



ym 



3Q-P , -3(R-S) 



3Q-P „-3(fl-S) 



2/i 



3Q+P _ ^ 3(R+S) ^ \ 



^(F + G+f + f) y 2 (F + G 



2 y 

„-3Z 

4 



^ ) 



-3ff± ^9? J 2 -8^J (3Af +4AJ ) 

There exists another nontrivial solution (x) = u x (l, 1, 1) with v x = 4(3A*+4A*) — 

solution is not of interest for our purposes. 



But this 
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All parameters appearing here are real, and equal to 

~ » 9 l + Ay 2 + y 2 2(1 - 2y 2 + yT) cos^i 1 + yl 

or a er 

2 2 , 1 + Ayf + yj (l-2y 2 + y 2 )cos^ 1 
F=l + y 1+ y 2 + — , G= , 

^ /i 2n cos ^i2 + 2a cos ^ 2 ^ ^ 2x sin ^i2 - Sasin^ 

Z, = (1 - !/2) ab ' Z =( 1 -f2) ^ , 

p = 1 + + yj _ + (1 - 2,; + rf) cos^ 
ar a 

o /-, 2n cos ^12 - a cos ^ 2 . 2 N sin ^i2 + «sin^ 2 
= ~ h > S = (l-y 2 ) , 

R = (l-2y 2 + y 2 )smjj 1 ^ 
In Eq. (160]) the parameters \l/2 are defined by 

= C^i + S ? 3 { ( F + G + X + ^) (C23 + ^ + ^ ^ + G ~ X) Sin ^ 23 } 

- y sin26 l i3|3c23(-R + 5 1 ) sin^cp - (3Q + P) cos^cp 
+ IJ2S23 (3(-R - 5) sin 5c*p + (30; - P) cos S C p) j 
* 2 = (P + G + — + — J (s 2 23 + y 2 2 c 2 23 ) -y 2 (F + G-—J sin 2fl 23 • (B2) 
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